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Abstract 

Massive spectral sum rules are derived for Dirac operators of SU(N C ) gauge theories with 
Nf flavors. The universal microscopic massive spectral densities of random matrix theory, where 
known, are all consistent with these sum rules. 
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The derivation of exact spectral sum rules for eigenvalues of the Dirac operator in QCD by Leutwyler 
and Smilga Q has recently led to a remarkable series of results. Most impressive are the developments 
from random matrix theory. There is now quite compelling evidence that the so-called microscopic 
spectral density of QCD (and hence, as a simple by-product also all spectral sum rules) can be 
computed exactly from an appropriate large- ./V matrix ensemble |2[ ||. The essential input is the 
existence of a non-vanishing spectral density p(X) at the origin A = 0, which, on account of the 
Banks-Casher relation^] p(0) = iVS/V , S = (^>ip), translates into the existence of a chiral condensate 
in QCD. From this perspective the microscopic spectral densities embody information about the 
symmetries of the Dirac operators alone, without any reference to the detailed dynamics of the gauge 
field interactions. If correct, this represents a major step towards the understanding and classification 
of spontaneous chiral symmetry breaking in QCD and related theories. In fact, the spectral sum rules 
of Leutwyler and Smilga are then seen to be of one very particular universality class, out of many. 
The field-theoretic generalization of the sum rules to other patterns of chiral symmetry breaking has 
been given in ref. Q. Furthermore, in a very impressive series of papers [| it has been shown 
how to translate the (conjectured) catalog of universality classes of random matrix models into gauge 
theory language (see, e.g., ref. || for a recent review). 

Evidence for the above scenario has been mounting after the initial observation that massless spectral 
sum rules of QCD could be reproduced by microscopic spectral densities computed in random matrix 
theory. First, it has been shown that the microscopic spectral densities of random matrix theory 
indeed do fall into distinct universality classes, independent of the detailed matrix model interactions 
0. Second, there are now direct numerical measurements of the microscopic spectral density for 
quenched SU{2) gauge theory, and the agreement with predictions from random matrix theory is 
quite spectacular Q. 

It was suggested by Shuryak and Verbaarschot || that the microscopic spectral density (and corre- 
sponding spectral sum rules), which normally are considered for massless Dirac operators, could be 
generalized to the massive case by an appropriate rescaling of masses. Such an extension is in fact 
essential for future comparisons with lattice gauge theory, where dynamical fermions in a finite vol- 
ume necessarily must be massive. Very recently, the associated double-microscopic spectral densities 
(called so because both eigenvalues and masses are rescaled as the volume N is taken to infinity) have 
been derived from random matrix theory, and proven to be universal as well ||]. A few examples 
showed agreement with corresponding massive sum rules of QCD. 

The purpose of this paper is to more systematically derive massive spectral sum rules for massive 
Dirac operators, and to confront these sum rules with the recent results []|] from random matrix 
theory. Throughout we restrict ourselves to the region 1/Aqcd iV 1 / 4 <C Here Aqcd is a 

typical hadronic scale in QCD, and is the pion mass. 

Gauge group SU(N C ),N C > 3, Nf fermions in the fundamental representation: 

The coset is here SU(Nf)L x SU(Nf)n/ SU(Nf), and the matrix model ensemble is that of chiral 
unitary matrices Let us start with the case Nf = 1, and the sector of zero topological charge: 
v = 0. Massive spectral sum rules can now be derived in several ways. Consider the general expres- 
sion for the partition function Z v in the topological sector v (I n (x) denotes the nth modified Bessel 

1 Here N denotes the volume of space-time. We use this notation for convenience, since it in matrix model language 
essentially corresponds to the size of the random matrices. 
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function) Q, 



H = NT,m. 
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From the expansion of the Bessel function, and a comparison with the corresponding formal expansion 
of the full QCD partition function it follows that 
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where the j v ,n s denote the zeros of the Bessel function J u (x), and where the zero mode has been 
omitted from the sum on the left hand side. 



The double-microscopic spectral densities p^ f \X; Hi,... ,PNf) are defined by 
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where p(A) is the ordinary, macroscopic, spectral density: 

P (A) = (£ 'S(X - X n )) 



Hi = rrtiNY, fixed 
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With the help of this definition the sum rule (y) can be written 
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We now compare this prediction with that of the chiral unitary random matrix ensemble. In ref. || 
the derivation of the double-microscopic spectral density was restricted to the sector of zero topological 
charge, v = 0. However, a simple analytical expression was also given for the case of one massive and 
Nf — 1 massless fermions, which conveniently can be written^ 
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Comparing with the general matrix model expression Q, this immediately gives us also the double- 
microscopic spectral density for Nf = 1 in the sector of arbitrary topological charge v: 
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Substituting the above expression into the sum rule (JSJ) , and performing the integral, one finds that 
this sum rule indeed is satisfied. 



One notices that the spectral sum rule ((B]) essentially is a rewriting of the (mass dependent) chiral 
condensate. The only slight complication is that the sum rules conventionally are written in terms 

2 Analogous expressions valid for any v for higher values of Nf are contained in the general solution given in ref. |?| 
by considering the case of v massless fermions and Nf massive ones. 



2 



on non-zero eigenvalues only. Because the partition function for one massive flavor can be written as 
(A n > and, for convenience in all that follows, v > 0) 



[dA„] v m v Y[(X 2 n + m 2 )exp{-S[A}} , 



(8) 



it follows that the general formula for the contribution of non-zero modes to the mass-dependent chiral 
condensate E(m) reads 

\*pzM-$ ■ (») 



In terms of the eigenvalues X n it becomes 
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This formula, and generalizations for higher inverse moments, will give us all the required massive 
spectral sum rulesj^j For example, we immediately recover the expression (|5|) without having to resort 
to the infinite summation over zeros of the Bessel function. It is also worthwhile noting that the 
massless spectral sum rules of Leutwyler and Smilga are recovered from the above expression in the 
massless limit: 
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We now turn to the general case of Nf > 2. The partition function Z v for equal masses m = ///(iVS) 
was found in the original work of ref. Q as the determinant of an Nf x Nf matrix, 



Z u = detM , 



1/ 
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More recently Jackson, §ener and Verbaarschot |^] have given the general formula for different masses: 
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is given in terms of derivatives of the Bessel function, and 

i<j 



is the Vandermonde determinant of masses. 

For equal masses, the generalization of eq. (|i~0| ) to an arbitrary number of flavors Nf is 
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The restriction to non-zero modes is not essential, and done here only in order to have a well-defined massless limit. 
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To illustrate this, consider the case Nf = 2, where this sum rule becomes 

/y / 1 \ iV 2 S 2 \ I u {^)I v+1 (y) - I u ^)I u+2 ^) \ 

\V >?n + ™ 2 l v ^ 1 W-I V+ l(jl)I„-l(j*) J ^ ^ 

As a simple by-product we recover the Leutwyler-Smilga sum rule for Nf = 2 in the massless limit [fij: 

The explicit expressions become increasingly involved with growing Nf. For Nf = 3, the sum rule 
analogous to IfFjf) reads 



/ 1 \ iV 2 S 2 

\ n A n+ m I v GflZ u 

-I v -\.{p)I v {p)I v+ 2{ji) + I u -i(fi) 2 - I v -2([j,)I u (n)I u+3 (n)) ,(19) 

where 

Z„ = ^(/j) 3 +^ + i(/j) 2 /^ 2 (^)+^-i(^) 2 /^ . (20) 

By explicitly taking the limit // — > on the right hand side of eq. (|l9|), we also here confirm that it 
correctly reduces to the massless sum rule ( |l8|) for Nf = 3. 

These massive sum rules can now be confronted with the predictions from random matrix theory. 
Taking degenerate masses, the double-microscopic spectral densities can be written, for u = 0: 

J2) f . lA _\C\\r (n 2, T (t) 2] 2|c|(M/iMJo(C) + qoMJi(C)) 2 (21) 

PS (C,M,M)- T [Jb(C) + M0 \ - (C 2 + /1 2)2 [/o(/l) 2_ /l(/l) 2 ] ( 21 ) 

for TVy = 2. Using this density to evaluate the left hand sides of eqs. ([l7|), one verifies that that this 
identity is satisfied. Even for equal masses, the expressions that must be integrated in order to verify 
the sum rules grow rapidly with increasing Nf. 

Many more different kinds of massive spectral sum rules can be derived from the partition function 
for different quark masses (|l3D . For example, for Nf = 2 and v = we get, from eq. (|l3|), 

/y , 1 \ = N 2 T? f /i(jm)/x 2 /i(/i2)-/o(/x 2 )m/o(/xi) | 2/xi 1 

\n A n + m l/ I , = 2/Xl I J (^1)^1(^2) - Io^2) f^lh(fJ-l) ix\-\l\]' 



We note that this correctly reduces to the sum rule (|17| ) in the limit of degenerate fermion masses 
(cancellations are delicate, and occur up to 2nd order). By numerical integration we have verified to 
high accuracy that the microscopic spectral density for two fermion flavors of different masses H , 



pS 2) (C;mi,m 2 ) = ^(j 2 (C) + Ji(0 : 



\cm-t4) [w/i(/xi)Jo(c)+ao(w)Ji(c)][M2A(/x2)^o(c)+ao(M2)Ji(c)} 



(C 2 +^i)(C 2 +Mi) /ii/i(/ii)/ (M2) - H2lo(ni)Ii(n 2 ) 

is consistent with this sum rule. 
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Massive spectral sum rules originating from higher derivatives of In Z v with respect to (different) 
fermion masses can be derived straightforwardly. Physically they are related to higher susceptibilities. 
Using the general expression for the universal massive double-microscopic spectral correlators of ref. 
p|, also these generalized sum rules are open to checks from random matrix theory. Higher massive 
sum rules for different fermion masses satisfy non-trivial identities as one or more of the masses are 
sent to infinity, due to decoupling. These identities have their direct counterparts in the consistency 
conditions for p$ and higher correlators that were explained in ref. ||. 

Gauge group SU(N C ),N C > 2,Nf fermions in the adjoint representation: 

The coset is SU(Nf)/ SO(Nf), and the relevant random matrix models belong belong to the chi- 
ral symplectic ensemble H, Massless spectral sum rules have been been derived by Leutwyler and 
Smilga for Nf = 1 and 2 JlJ , and later generalized to arbitrary Nf by Smilga and Verbaarschot Q : 

N 2 T? , . 

(24) 




4(* + (iV> + l)/2) • 

Here v = N c v is the number of zero modes, and on the left hand side the sum runs over the doubly- 
degenerate, positive eigenvalues, counted once. The partition function has been given explicitly for 
Nf = 1 yj, where Z u = Ip(n). The massive spectral sum rules therefore coincide with those of 
SU(N C > 3) for Nf = 1 and v replaced by v. In particular, eq. (||), with v replaced by D, for the 
simplest of these massive sum rules. Also the partition function for Nf = 2 and degenerate masses is 
known jjj]: 

Zv ~ ^ n\(n + 2u)\(2n + 2U + 1) ' ^ 25 ^ 
After some manipulations, we have simplified this expression to 

Z = 7 (2 M ) + | [/ (2/i)L 1 (2/x) - / 1 (2/x)L (2/x)] (26) 

for P = 0, and 

Zu = (-If {Wfi) + - [J (2m)L 1 (2 M ) - Ii(2/i)L (2 M )] - - ^(-l) fc / 2fe+1 (2^)) (27) 

for v > 1. Here L n (x) is the nth modified Struve function. We note that the right hand side of eq. 
(p7|), despite appearances, is positive also for odd v. 

Massive spectral sum rules can now be derived by straigthforward differentiation, but for large values 
of v they yield rather unwieldy expressions. We restrict ourselves here to displaying only the simplest 
sum rule for v = 0, which in fact simplifies considerably: 

Using the expansions Lo(x) = 2x/n + . . . and Li(x) = 2x 2 /(3tt) + . . ., one verifies that this massive 
spectral sum rule reduces to the massless spectral sum rule ( p4| ) in the limit //— >0. There are presently 
no predictions from random matrix theory with which to compare the massive spectral sum rule 
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Staggered fermions of SU(2) lattice gauge theory have the symmetries of the chiral symplectic matrix 
ensemble ||. The massive sum rules for Nj = l and Nf = 2 shown above may therefore soon be tested 
by Monte Carlo simulations. 

Gauge group SU(2),Nf fermions in the fundamental representation: 

The coset is here SU(2Nf)/Sp(2Nf), and the matrix models are of the chiral orthogonal ensem- 
ble ||. The massless spectral sum rules have been derived by Smilga and Verbaarschot 0|, and found, 
for the simplest case, to be of the form 

\\ K/ v " 4(^ + 2^-1) • (29) 

The partition function Z v in the sector of topological charge v can, for Nf flavors of equal masses 
m = fj,/(N'E), be written in terms of a Pfaffian 

Z» = (2 Ar-l)!! P/(A) ' Aij = 0'-0^+i+i(A*) • (30) 



For Nf = l this formula gives Z v = and the massive spectral sum rules therefore coincide with 

those of SU(Nf > 3), JVy = 1, and in particular eq. (g) for the simplest sum rule. For Nf = 2 we get, 
from eqs. (|l6|) and (|30|), 



1 \ iV 2 S 2 f 2I u ^)I u+1 (fi) - 3J„-i(//)J„ +2 0z) + J^aQ^+aQi) ■ 
,V A n + W„ ^ I 3J |/ (M) 2 -4/ I/ _ 1 (^)J H . 1 0*) + / I/ _ 2 (/i)J H . 2 ( M ) I' 



As a direct check of this expression, we also here find that it reduces to the known result (|29| ) in the 
massless limit fi — > 0. Expressions for larger values of Nf can be worked out analogously. There are 
presently no predictions from random matrix theory with which to compare these massive spectral 
sum rules. 



We end this paper with a few comments on the physical significance of these massive sum rules. The 
whole idea is to rescale both eigenvalues A and fermion masses mj at the same rate as the volume is 
taken to infinity, so that £ = AiVE and \ii = rriiNT, are kept fixed. In addition, we are restricted to 
volumes N that satisfy 1/Aqcd "C iV 1 / 4 <^ l/m n . To obtain a chiral condensate of the theory with 
massless fermions in the infinite- volume limit, it is well-known that one should first take the volume 
N to infinity, and then take m; to zero. This is not the limit considered here, where the m; — > 
at precisely just the proper rate with N — > oo to leave a well-defined double-microscopic massive 
spectral density. Highly non-trivial information about the finite-volume correlations of eigenvalues is 
contained in this limiting function. From the point of view of lattice gauge theory simulations, the 
limit we have discussed here is completely realizable. The required tuning of quark masses as different 
lattice volumes are considered is an entirely natural and feasible operation in Monte Carlo simulations. 
There are therefore good reasons to believe that the massive spectral sum rules discussed here may 
be tested in lattice gauge theory. It is already a highly non-trivial fact that the double-microscopic 
massive spectral densities of random matrix theory, in those cases where known, are consistent with 
these sum rules. 
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